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Abstract 

Vapour-liquid equilibria (VLE) and the influence of an inert carrier gas on homogeneous vapour 
to liquid nucleation are investigated by molecular simulation for quaternary mixtures of carbon 
dioxide, nitrogen, oxygen, and argon. Canonical ensemble molecular dynamics simulation using 
the Yasuoka-Matsumoto method is applied to nucleation in supersaturated vapours that contain 
more carbon dioxide than in the saturated state at the dew line. Established molecular models are 
employed that are known to accurately reproduce the VLE of the pure fluids as well as their binary 
and ternary mixtures. On the basis of these models, also the quaternary VLE properties of the bulk 
fluid are determined with the Grand Equilibrium method. 

Simulation results for the carrier gas influence on the nucleation rate are compared with the 
classical nucleation theory (CNT) considering the 'pressure eff"ect' [Phys. Rev. Lett. 101: 125703 
(2008)]. It is found that the presence of air as a carrier gas decreases the nucleation rate only 
slightly and, in particular, to a significantly lower extent than predicted by CNT. The nucleation 
rate of carbon dioxide is generally underestimated by CNT, leading to a deviation between one and 
two orders of magnitude for pure carbon dioxide in the vicinity of the spinodal line and up to three 
orders of magnitude in presence of air as a carrier gas. Furthermore, CNT predicts a temperature 
dependence of the nucleation rate in the spinodal limit, which cannot be confirmed by molecular 
simulation. 
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1. Introduction 



If significant reductions in greenhouse gas emissions are to be achieved, it is essential to min- 
imize the energy required for separating carbon dioxide from air as a prerequisite for rendering 
purification, r ecovery and sequestration processes economically and ecologically more efficient 
(|KauleL |2002|) . This requires an accurate understanding of the phase behaviour of carbon dioxide 
(CO2) in the earth's atmosphere under diverse conditions, corresponding to the respective technical 
applications. 

Condensation processes, initiated by nucleation, as well as vapour-liquid equilibria (VLE) are 
particularly relevant in this context, since they characterize capture and storage of CO2, i.e. phase 
separation and phase equilibrium. The present study concentrates on understanding homogeneous 
nucleation and VLE properties for quaternary mixtures consisting of CO2, nitrogen (N2), oxygen 
(O2), and argon (Ar) by means of molecular simulation. Thereby, the properties of the investig- 
ated fluid itself can be fully isolated from phenomena induced by impurities or boundary effects 
in the vicinity of a solid wall, which would be much harder to accomplish in an experimental 
arrangement. 

Although a systematic discussion of CO? nucleation was, surpri singly enough, published for 
the ambient conditions prevailing on Mars (^Ma atanen et al. . 2005 ). to the authors' knowledge 
no analogous study is available for the ecologically and technically more relevant atmosphere 
composition of our own planet. The present work closes this gap and characterizes the air pressure 
effect on the condensation process in a vapour containing more CO2 than at saturation on the basis 
of molecular models that are known to accurately reproduce VLE properties. 

Direct MD simulation of nucleation in the canonical ensemble with the Yasuoka-Matsumoto 



(YM) method is an established approach (lYasuoka and Matsumotolll998 . Matsubara et al.Ll2007|) 



It can be successfully applied to the regime where the supersaturation is sufficiently high to permit 
significant droplet formation in a nanoscopic volume within a few nanoseconds, but not as high 
that the nucleation rate /" is affected by deplet ion of the vapour du e to droplet formation within 
the same time interval ( Horsch et al. . 2008aubl . Chkonia et al. . 2009*). Lower nucleation rates can 



be de termined by forward flux sampling or methods based on umbrella sampling ( Valeriani et al. 



20051), whereas in the immediate vicinity of the spinodal line, where the metastable state would 
otherwise break down, a stationary value fo r the nucleation rate can be obtained by grand canonical 
MD simulation with McDonald's daemon (IHorsch and Vrabed . l2009h . 

The present work focuses on the intermediate regime, where the YM method is viable, and 
compares simulation results with theoretical predictions on the basis of two variants of the classical 
nucleation theory (CNT). 



2. The pressure effect 

If the pressur e effect (PE) is taken i nto account, the free energy G of formation according to 
CNT is given by dWedekind et al.L l2008h 



dGpE = lyUs - // - ^1—^ I dN + ydF, 



(1) 
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for a droplet containing N molecules, whereas the more commonly used - but thermodynamically 
inconsistent - simplified classical (SC) variant of CNT implies 



JGsc = (jUs-M) dN + ydF. 



(2) 



Therein, yUj and as well as /i and p represent the saturated and supersaturated values of chemical 
potential and pressure, respectively. While refers to the pure substance vapour pressure of the 
nucleating component, p is the total pressure including the pressure contribution of an inert carrier 
gas, if present. The droplet surface tension is given by y, the area of the surface of tension is F, and 
p' represents the density of the liquid which is assumed to be incompressible. Where the chemical 
potential difference Au = jx- jx^ appears in the SC expression, PE applies an 'effective' difference 



(IWedekind etal.Ll2008h 



p' 



(3) 



which is always smaller than Ajx. Assuming a spherical shape for the emerging droplets, i.e. 



dF = 



32n 



1/3 



dN, 



(4) 



and neglecting the curvature dependence of the surface tension, the free energy barrier of the 
nucleation process according to the PE variant of CNT 



2/3 



To, 



(5) 



is reached for the critical droplet size 



32n 
3^ 



Jo 

Ap^ 



(6) 



These expressions use the surfac e tension yn of the planar vapour-liquid interface as well as the 'in- 
ternally consistent' approach of iBlander and Kata (119721) which equates single-molecule droplets 
with vapour monomers and assigns them a free energy of formation AG = 0. The same results 
follow for the SC variant, with Ap instead of Ap^ in Eqs. ([5]) and Q. 

The long-term growth probability (3(M) of a droplet containing M molecules can be given by 



Q{M) 



Jl 



/2AG(^ 



-1 



(7) 



under the approximation that the r eaction coordinate o f the nucleation process only depends on 
the droplet size order parameter N (|Horsch et al.Ll2009r) . The transition rate of vapour monomers 
of the nucleating component through an interface, normalized by the surface area of the interface, 
is 

j3 = 'Pi2mnakTr^'\ (8) 
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according to kinetic gas theory, where k is the Boltzmann constant, nio is the molecular mass and 
P the pressure of the nucleating component in its p ure gaseous state a t the same partial density. 



The isothermal nucleation rate according to CNT is (|Feder et all Il966r) 



jr = ZF*l3pexp 



-AG* 
kT 



(9) 



which depends on the surface area F* of the critical droplet, the lZel'dovid (|l942h factor 



Z = 



InniokT 



(10) 



N=N* 



as well as the density of vapour monomers belonging to the nucleating component p. Droplet 
overheating due to rapid growth, however, is neglected in the expression for jj- The monomer 
density in the metastable state can be obtained from 



(11) 



a series that usually converges very fast, wherein p is the total density of the supersaturated vapour 
and Yq is the mole fraction of the nucleating component in the supersaturated vapour. The presence 
of a carrier gas also influences the thermalization of growing droplets, facilitating the heat transfer 
from the liquid to the surrounding vapour and thereby decreasi ng the amount of ov erheating. This 



efi'ect is covered by the thermal non-accomodation prefactor of iFeder et al.l (|1966|) 



&2 



+ q 



:JT, 



consisting of 



(12) 



(13) 



N=N* 



as well as the mean square fl uctuation of the kinetic energy for the vapour molecules (IFeder et al 



1966LIWedekind et al.Ll2008h 



kT^ 



k\^ Yim'l\c,j + kl2) 



VUYom]'\c,,o + kl2) 



(14) 



where h' and h" are saturated liquid and vapour enthalpy, respectively, c,,,, is the isochoric heat 
capacity of component i, where / = indicates the nucleating component and I < i < K the 
components of the inert carrier gas. Within the scope of the present study, the heat capacity of 
the pure saturated vapour is used for Cv,o whereas for the other c,,,, the value in the limit of infinite 
dilution is used, since CNT assumes the carrier gas to have ideal properties (tWedekind et al.L 
2008h . 
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p' \kTl 



(15) 



Within the framework of CNT, it follows that the prefactor 

does not depend on the pressure contribution of the carrier gas, and neither does /3, while the 
influence of the carrier gas pressure on p, which is usually similar in magnitude to pY^, is of minor 
significance. Th is eliminates all contributions to the pressure effect except for those discussed by 
Wedekind et aP (E0O8,) as 

JPE ''pE JT,PE 



Jr,sc t'pE + <?PE Jt,sc 



(16) 



Norma lized to unity for the pure fluid (Yq 
(|2008|) pressure efi'ect can be expressed as 



1) with the PE variant of CNT, the lWedekind et al 



"WiYo) 



^pe(^o) [^^e(^o = 1) + ^pe(>^0 



[bl^{Yo) + ql^(Yo)] bl^(Yo = l) 



1)] ^ _^/ AG*^(y o = l)-AG 



PE 



kT 



(17) 



given that 7V,sc> the denominator of Eq. (fT6l) . does not depend on Fo- 
under certain cond itions, the p r essure efl'ect does not exceed the experimental uncertainty and 



can thus be neglected (|Iland et al.L |2004|) . In other cases, however, the influence can be experi- 



mentally detected, with apparently contradictory results: sometimes / inc rease s with the amount 
of car rier gas, in other cases the opposite tendency is observed (Hyvarinen et al.Ll2006llBrus et all 



20081). The 'W factor explains this in principle, since it combines the thermal non-accomodation 
factor, which increases with Yq 0, and the free energy effect that leads to an effective chem- 
ical potential difference Afi^ < Ai x. Depending on the th ermodynamic conditions, each of these 



contributions can be predominant (IWedekind et al.Ll2008|) . 



The main inaccuracies of CNT concern the surface tension as well as the droplet surface area. 
For the surface tension, deviations from t he capillarity approximatio n y ~ jo are known to oc- 
cur for nanqscopica lly cu rved interfaces (|Moody and Attardl boM ISzybisz and Urrutial 120031 



Horsch et al.L l2008b|) . The iTolmanI (|1949|) approach implies huge deviations from Jq for droplets 
on the molecular length scale corresponding to the Tolman length 6. In particular, for > 6 the 
dependence of y on the surface of tension radius "Ry can be expressed as 



7 = 70 exp 



(18) 



whereas y ~ Ky becomes valid for Ky «; 6. From the Laplace equation along with the iTolman 



approach it can be deduced that the area of the surface of tension is given by 

2dN 



dF = 



for an incompressible fluid, where the relation 



'Ry — 



3N 
Anp' 



1/3 



5, 



(19) 



(20) 



follows for the dependence of the surface of tension radius on N. Overall, this leads to a signific- 
antly increased surface area for droplets on the molecular length scale. 
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3. Simulation methods and models 



The evaluation of the theoretical predictions relies on knowledge about the chemical potential 
difference between the saturated and the supersaturated state. This was obtained for pure CO2 by 
Gibbs-Duhem integration 

AiuiT,p)= r (21) 
Jps(r) p(T,p) 

using MD simulation results of small systems (A^ x 10 000) in the metastable vapour regime. The 
carrier gas influence according to the presented variants of CNT was evaluated by assuming ideal 
gas properties for air as well as ideal mixing behaviour, i.e. 



p = {l-Y,)pkT + 'P. 



(22) 



For the homogeneous nucleation simulations, the YM method was applied to relatively large, 
but still nanoscopic systems with NiCOj) = 300 000. The total number of molecules was up 
to = 900 000 such that the carrier gas with A^(N2) : A^(02) : N(Ar) = 7812 : 2095 : 93 
corresponded to the earth's atmosphere composition. The condensation process is thereby re- 
garded as a succession of three char acteristic stages: relaxation, nucleation, and droplet growth 
(|Chesnokov and Krasnoperovl |2007|) . During the nucleation stage, the droplet formation rate 
I(M), i.e. the nur nber of droplets containing at le ast M molecules formed over time, is approx- 
imately constant (lYasuoka and Matsumotol Il998h . The droplet formation rate depends on the 
threshold size M and is related to the nucleation rate by 



I(M) 
V 



J 



Q{M) 



(23) 



since Q{M) indicates the probability for a droplet containin g M molec i iles to reach macroscopic 
size. Liquid and vapour were distinguished according to a IStillingeii (1 19631) criterion such that 
molecules separated by distances of their centres of mass below 5.08 A were considered as part 
of the liquid. Biconnected components, where any single connection can be eliminated without 
disrupting the internal connectivity, were defined to be liquid droplets. 

Molecular models for Ar, which can be represented by one Lennard- Jones (LJ) site, and for 
CO2, N2 as well as O2, which can be represented by two LJ sites separated by the elongation £ with 
a superimposed quadru pole moment O in th e molecule's centre of mass (2CLJQ), were adjusted to 
pure fluid VLB data by lVrabec et al.l(|2001r) . cf. Tab.[T] If adequate values f or the unlike d i spersiv e 
interaction energy are used, so that binary VLB are reproduced correctly ('Vrabec et al.'. ' 2009ah 



cf. Fi g.[Tl ternary mixtures are accurately described without any further adjustment (Huang et al. 



2009h . In Tab. [T] the unlike energy param eters are indicate d in terms of the binar y interaction 
parameter ^ of the modified Bertheloti ( 1898 ) combining rule ( Schnabel et al. . 2007b ) 



^AB - ^(Sa^b) 



1/2 



(24) 



while the unlike LJ size parameter is determined as an arithmetic mean according to the iLorentz 
(Il881|) combining rule. This approach has also been vali dated with an emphas is on CO2 in particu- 
lar, confirming its viability for mixture s with N2 and O2 (| Vrabec et al.l . l2009b|) as well as hydrogen 



bonding fluids (|Schnabel et al.L 12007 ah . 





o 

type (T [A] 


e [meV] 


^A] 


Q [eA2] 


^(Ar) ^(02) ^(N2) 


C02 


2CLJQ 2.9847 


11.394 


2.418 


0.78985 


0.999 0.979 1.041 


N2 


2CLJQ 3.3211 


3.0072 


1.046 


0.29974 


1.008 1.007 


02 


2CLJQ 3.1062 


3.7212 


0.9699 


0.16824 


0.988 


Ar 


LJ 3.3967 


10.087 









Table 1: Molecular model para meters of IVrabec et al.l (1200 Ih and binary interaction parameters ^ adjusted to binary 
VLE data dVrabec et al.Ll2009ah . 



On that basis, quaterna ry phase equilibria were determined using the Gr and Equilibriu r n metho d 
(|Vrabec and Hassd . |2002|) . introducing Ar into the the system studied by IVrabec et al.l (|2009b|) . 
The Grand Equilibrium method calculates the vapour pressure as well as all dew line mole frac- 
tions ji by simulation for a specified temperature T and a specified bubble line mole fraction Xi 
for all components of the mixture. Although no experimental data are available for the quaternary 
mixture, the simulation results can be trusted due to the extensive validation of the models with 



respect to the VLE behaviour for all of the six binary (IVrabec et al.L l2009a|) and two of the four 
ternary subsystems (|Huang et al.ll2009|) . i.e. N2 + O2 + Ar as well as CO2 + N2 + O2. 
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Figure 1 : (Colour o n the web, b/w in print.) Experimental data (★) of jDodgg (1 19271) and simulation results (o) of 
Vrabec et al.l ( l2009ai) . using the Grand Equilibrium method with the molecular models given in Tab. [Tj for VLE of 
binary mixtures containing nitrogen and oxygen at temperatures of T = -193, -168, and -153 °C. 



4. Simulation results 

Grand Equilibrium simulations of the quaternary mixture CO2 -1- N2 -1- O2 -i- Ar were conducted 
for VLE covering a broad temperature range with CO2 bubble line mole fractions x{C02) of 0.910, 
0.941, and 0.969, cf. Tab.[2l Except for the highest temperature, which corresponds to 93 % of the 
critical temperature Tc for pure CO2, the mole fractions j(N2), 3^(02), and y(Ar) on the dew line are 
one order of magnitude higher than the corresponding bubble line mole fractions. This confirms 
that for temperatures sufficiently below rc(C02), air only accumulates to a limited extent in the 
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liquid phase. As a first approximation, it can therefore be treated as a carrier gas for CO2 nucleation 
so that the PE variant of CNT with a single nucleatin g component can be appl ied, as opposed 
to more complex mixtures such as ethanol + hexanol dStrey and ViisanenL Il993|) . water-alcohol 



mixtures (tViisanen et al.L 11994 IStrey et al.l . 119951) . or water + nonane + butanol (INellas et al. 



20071), where multi-component nucleation occurs. 

From MD simulation of small metastable systems, the spinodal value 5 ^ of the supersaturation 
with respect to density, which is defined as 



S = 



p"(Ty 



(25) 



wherein p"{T) is the saturated vapour density of pure CO2, was determined to be in the range 
4.3 < 5^ < 5.1 at -44.8 °C, 3.6 < 5^ < 4.3 at -34.8°C, and 3.0 < 5^ < 3.6 at -23 °C for = 1. 
At these temperatures, canonical ensemble MD simulations for CO2 nucleation were conducted 
using the YM method with CO2 mole fractions of 1/3, 1/2, and 1 at supersaturations below the 
spinodal value S^, but still high enough to obtain statistically reliable droplet formation rates in a 
nanoscopic volume on the timescale of a few nanoseconds. 



YM droplet formation rates / from the present work as well as a previous study (|Horsch et al.L 
2008 a|) are shown in Fig. [2] for 70 = 1^ i-e- pure CO2. The dependence of / on the threshold size 



M reproduces the typical picture: for low threshold sizes (probably smaller than N*) the droplet 
formation rate can be elevated by several orders of magnitude, and it converges for M » N* 
under the condition that the depletion of the vapour can be neglected (|Yasuoka and Matsumoto . 
19981). For very large values of M - not shown in Fig. |2]- the droplet formation rate decreases 
again, because the presence of many large droplets implies that a substantial amount of the vapour 
monomers have already been consumed by the emerging liquid phase. 

The most striking observation is that while both variants of CNT predict the value of j in the 
spinodal limit to increase with temperature - mainly because T occurs in the denominator of the 
exponential in the Arrhenius term of Eq. ^ - the simulation results do not exhibit any significant 
temperature dependence for the attainable value of j. In the spinodal limit, the nucleation rate 
appears to be about j(r, 5^) « 10^^ cm^^s ' over the whole temperature range. The pressure 
eff'ect in the pure fluid, expressed by jsd Jpe, is most significant at high temperatures, because this 



T[°C] x(C02) 


pJMPa] y/x(N2) y/x{02) y/x(Ar) p' [mol/1] p" [mol/1] h" - h' [kymol] 


-90.3 0.969 
0.941 
0.910 


2.53(8) 40(2) 24.5(9) 24.7(9) 29.24(1) 1.854(2) 17.092(9) 
3.9(1) 21.5(8) 13.9(5) 13.7(5) 29.13(1) 3.031(5) 16.27(1) 
6.0(2) 13.7(4) 9.1(3) 9.0(3) 29.01(2) 5.29(2) 15.12(1) 


-40.3 0.941 
0.910 


4.38(5) 14.0(2) 10.8(1) 10.3(1) 24.94(3) 2.663(5) 13.30(1) 
5.30(4) 9.7(1) 7.72(8) 7.34(8) 24.47(2) 3.247(7) 12.64(1) 


9.7 0.969 
0.941 


5.97(4) 5.6(1) 4.88(8) 4.72(7) 19.3(2) 3.90(1) 8.55(5) 
6.98(3) 4.31(5) 3.82(5) 3.67(4) 18.55(9) 4.63(2) 7.79(4) 



Table 2: VLE data for the quaternary system CO2 + N2 + O2 + Ar. The Hquid composition is equimolar in nitrogen, 
oxygen, and argon, i.e. x(N2) = x{02) = Jt:(Ar) = [1 - x(C02)] / 3, and values in parentheses indicate the uncertainty 
in terms of the last given digit. 
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corresponds to a lower density of the liquid and because Afx is smaller so that the relative deviation 
between A/u^ and A/u is increased. 
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Figure 2: (Colour on the web and in print.) Pure CO2 nucleation rate j according to the PE ( — ) and SC ( ) 

variants of CNT in comparison to IjV for threshold sizes of M = (A) 50, (o) 75, and 250 (•) molecules from canonical 
ensemble MD simulation over supersaturation S - p/p"{T) at temperatures of 7" = -44.8, -34.8, -23, and -4.2 °C. 
The simulation results for T - -4.2 °C are taken from previous work and were obtained using a different cluster 



criterion ( Horsch et al. , 2008a ) 



Table [3] indicates the results for the carrier gas effect on CO2 nucleation at T = -44.8 °C. As 
usual, / decreases when larger values of the threshold size M are regarded, but it can also be seen 
that this effect is clearly stronger when more air is present in the system. This leads to values for 
the overall carrier gas effect ^ on the droplet formation rate 



I(S,Yo,M) 
7(5,70 = 1,M)' 



(26) 



that are greater than unity for relatively small values of M, but converge to values significantly 
below unity as M is increased. This result can be understood if the carrier gas effect on N* ac- 
cording to Eqs. ([3]) and ^ is considered: with a higher total pressure, Afx^ decreases which affects 
the critical droplet size to the third power, leading to significantly larger values of N^^, cf. Tab. |3l 
Thus, for relatively small threshold sizes, the long-term growth probability is significantly smaller, 
cf. Eq. (|7]), which in turn increases the droplet formation rate according to Eq. (|23]) . Hence, the ap- 
parently contradictory values of ^ are actually consistent and correspond to a negative dependence 
of the nucleation rate on the carrier gas density. 

While this qualitatively confirms CNT with the pressure effect, which leads to 'W factors 
on the order of 0.1 for Yo = 1/2 and 0.01 for Yq = 1/3, the actual decrease of I approaches 
the range 0.3 to 0.4 in the limit of large threshold sizes for both values of Yq. This impression 
consolidates itself if the results for T = -34.8 and -23 °C are also regarded, cf. Fig.[3]and Tab. lU 



The normalized I Wedekind et al.l (|2008|) pressure effect 'W, corresponding to the deviation between 
7PE(io) and 7pe(5^o = 1)> decreases even faster with Yq —> at high temperatures. Qualitatively, 
this is confirmed by simulation results, e.g. no nucleation at all was detected at these temperatures 
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Table 3: Droplet formation rate and car rier gas effect from YM canonical ensemble MD simulation as well as critical 
droplet size (in molecules), normalized Wedekind et alJ (|2008) pressure effect, nucleation rate and droplet formation 
rate according to the PE variant of CNT, in dependence of supers aturation and mole fraction of CO2 in the vapour as 
well as the YM threshold size (in molecules) at a temperature of T = -44.8 °C. 



for Fq = 1 /3, which implies j < 10^^ cm ^ However, from the available results for = 1/2 it 
is evident that the pressure effect is overestimated by CNT, in particular atT = -23 °C. 



5. Conclusion 

Vapour-liquid coexistence in equilibrium and non-equilibrium was studied by molecular simu- 
lation for systems consisting of CO2, N2, O2, and Ar. For the nucleation of pure CO2, it was found 
that both the SC and the PE variant of CNT underestimate the nucleation rate by up to a factor 20 
for SC and between one and three orders of magnitude for PE. 

It should be noted that this result for pure CO2 nucleation is both qualitatively and quantitat- 
ively similar to the deviation of CNT for homogeneous nucleation of the truncated-shifted LJ fluid 
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Figure 3: (Colour on the web and in print.) Nucleation rate y'pE according to CNT with the pressure effect in 

comparison to I jV for threshold sizes of M = 50 (A ■ ■ ■ A), 100 (□ □), and 150 (■ — ■) molecules from 

canonical ensemble MD simulation in dependence of Y^^ - 1, i.e. the ratio between air and CO2 molecules in the 
system, for T = -44.8 °C with S = 3.42 and T = -34.8 °C with S = 2.8 as well as T = -23 °C with S = 2.34. Note 
that in the latter two cases, no nucleation was detected during the MD simulation run for ' - 1 = 1 and 2, implying 
that//y < 10^5 cm-^s-^ 




(jHorsch et al.L 1200 8b|) . In that case, it was shown that the deviation can be corrected by a thermo- 
dynamically consistent approach that takes both the curvature dependence of y and the increased 
surface area due to fluctuations of the droplet geometry into account. Future work could elaborate 
on that and validate whether a theory based on these considerations, e.g. as expressed by Eqs. (fTSl) 
and (fT9l ). applies for CO2 as well. 

The effect of air as a carrier gas on CO? nucleation as deter mined from MD simulation qualit- 
atively confirms the theory outlined by IWedekind et al.l (|2008|) . Although significant quantitative 
deviations were found, these are partly due to the fact that at the high densities corresponding to 
the present simulations, air cannot be reliably described by the ideal gas equation. The non-ideality 
leads to a lower total pressure and thereby reduces the magnitude of the pressure efl'ect to a certain 
extent. 
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Table 4: Droplet formation rate and carrier gas effect from YM canonical ensemble MD simulation as well as critical 
droplet size (in molecules), normalized Wedekind et al. (2008) pressure effect, nucleation rate and droplet formation 
rate according to the PE variant of CNT, in dependence of temperature, supersaturation and mole fraction of CO2 in 
the vapour as well as the YM threshold size (in molecules). 
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